Three-dimensional two-pion source image from Pb+Pb collisions at y/sNN=17 .3 GeV: 

new constraints for source breakup dynamics 
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Source imaging methodology is used to provide a three-dimensional two-pion source function 
for mid-rapidity pion pairs with px < 70 MeV/c in central (0 — 7%) Pb-fPb collisions at 
-ysjy^=17.3 GeV. Prominent non-Gaussian tails are observed in the pion pair transverse momentum 
(outward) and in the beam (longitudinal) directions. Model calculations reproduce them with the 
assumption of Bjorken longitudinal boost invariance and transverse flow blast-wave dynamics; they 
also yield a proper time for breakup and emission duration for the pion source. 

PACS numbers: PACS 25.75.Ld 



Over the last several decades, the expansion dynam- 
ics and breakup lifetime for the exotic matter produced 
in relativistic heavy ion (RHI) collisions, have been of 
paramount interest [l|, Q. Such enormous energy den- 



sities are created in the RHI collision zone, that decon- 
finement of nuclear matter is expected To gain an 
understanding of this state of matter, it is essential to 
study its dynamical evolution. The space-time extent of 
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hot particle emission sources in heavy ion colhsions has 
been studied for years via final-state correlations between 
two particles [1]. 

Years ago, pioneering work at the Alternating Gradi- 
ent Synchrotron (AGS) ^ and by the NA49 Collabora- 
tion at the CERN Super Proton Synchrotron (SPS) 0, 
exploited the Hanbury-Brown Twiss (HBT) correlations 
of hadron pairs in conjunction with fits to particle spec- 
tra, to estimate the dynamical properties of the reaction 
source in a blast wave model. The NA49 data for cen- 



Detector achieves precision tracking and particle iden- 
tification using four large Time Projection Chambers 
(TPCs). The first two of them are mounted in precisely 
mapped magnetic fields with total bending power of up 
to 9 Tm. Charged particles are detected by the tracks left 
in the TPCs and are identified by the energy deposited in 
the TPC gas. Mid-rapidity particle identification is fur- 
ther enhanced by a time-of- flight wall (resolution 60 ps). 
Charged particle momenta are determined from the de- 
flection in the magnetic fleld. With the NA49 setup. 



tral Pb-|-Pb collisions at y/s 



NN 



17.3 GeV indicated a a resolution of 6p/p^ w (0.3 - 7) x 10"* (GeV/c) 



strong longitudinal flow with an approximately boost in- 
variant longitudinal velocity profile. These data also sug- 
gested a transverse expansion of the pion emission source 
with a duration of 8-9 fm/c. A more recent analysis 
confirms these earlier findings while extending the beam 
energy dependence of the measurements to five separate 
SPS energies. 

Several years ago, an alternative technique based on 
source imaging was introduced for model-independent 
extraction of emission sources in the pair-center-of-mass 
system (PCMS) [1, S This new methodology has 
provided a more faithful and detailed extraction of the 
actual ID source function [Til [l^ . Recent theoretical 
developments [l^, [13, 0, Il6j enable the extraction of 
three-dimensional (3D) profiles of the emission source. 

This methodology, in both its ID and 3D forms, 
has been employed for Au-|-Au reactions at v^TVAf ~ 
200 GeV to obtain detailed information on both the spa- 
tial and the lifetime extents of the created emission source 
[l3]- Here, we use the new methodology to again 
study reaction dynamics at the SPS but with identified 
pion correlations and extensively developed imaging tech- 
niques that explicitly include Coulomb effects and do not 
assume Gaussian sources. The resulting non-Gaussian 
source functions are interpreted in the context of a, power- 



ful new simulation model, THERMINATOR 
This approach explicitly includes all known resonance de- 
cays, longitudinal expansion, transverse expansion and a 
freeze-out hypersurface. 

In this study, the source imaging technique is used to 
analyze NA49 Collaboration data for central (0 — 7%) 
Pb+Pb collisions at y/sj^j^ = 17.3 GeV, obtained at the 
SPS. Model comparisons allow tests of different aspects 
of the dynamics and, in particular, the extraction of the 
proper time for breakup and emission duration for the 
pion emission source. The picture that emerges has many 
similarities to that from the early work but also 

adds interesting features and conclusions that preclude 
direct comparison. 

The data presented here were taken by the NA49 Col- 
laboration during the years 1996-2000. Lead beams of 
158^ GeV from the CERN SPS accelerator were made 
incident on a lead foil of thickness 224 mg/cm^. De- 
tails of the experimental setup are discussed in Refs.[7] 
and |2lj. Briefly, the NA49 Large Acceptance Hadron 



achieved. Event centrality is determined by a forward 
veto calorimeter which measures the energy of spectator 
matter. Approximately 3.9 million central events were 
analyzed. 

The 3D correlation function, C(q), and its ID angle- 
averaged counterpart C(g), were obtained as the ratio 
of pair to uncorrelated reference distributions in rela- 
tive momentum q, for n~^Tr~^ and 7r~7r~pairs. Here, 
q = (pi~P2) jg ija,ff of the momentum difference between 
the two particles in the PCMS, and q is the modulus 
of the vector q. The pair distribution was obtained by 
pairing particles from the same event; the uncorrelated 
distribution was obtained by pairing particles from dif- 
ferent events. The Lorentz transformation of q from the 
laboratory frame to the PCMS is made by a transfor- 
mation to the longitudinally co-moving system (LCMS) 
frame along the beam direction followed by a transforma- 
tion to the PCMS along the pair transverse momentum 
[13. C{q) is observed to be flat for 50 < g < 100 MeV/c 
and is normalized to unity over this range. 

Mid-rapidity {\yL — yo\ <0.35, where ul and yo are par- 
ticle and nucleus-nucleus centre-of-mass rapidities in the 
laboratory frame), low fc^ {kr < 70 MeV/c, where fcr is 
half the transverse component of the pair total momen- 
tum) TT+TT^ and TT^TT^ piou pairs were selected for this 
study. Track merging and splitting effects were removed 
by appropriate cuts on both the pair and uncorrelated 
distributions. The merging cut required the two parti- 
cles in a pair to be separated by at least 2.2 cm over 50 
pad rows in the vertex TPCs [7|. A 20% increase in this 
minimum separation has resulted only in changes within 
the statistical errors. Similar evaluations for the other 
cuts indicate an overall systematic uncertainty which is 
comparable to or smaller than the statistical uncertainty. 

The effects of track momentum resolution were as- 
sessed by jittering the momenta of the tracks in the data 
by the maximum momentum resolution, Sp/p^ « 7x 10^** 
(GeV/c)~^. The resulting re-computed ID and 3D cor- 
relation functions, which incorporate twice the effect of 
the momentum resolution, were found to be consistent 
with those obtained without momentum smearing. The 
correlation functions without additional smearing serve 
in the following as a basis for the extraction of source 
functions via imaging and fitting. 

The imaging procedure employed uses the ID imaging 
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code of Brown and Danielewicz 0, , which has been 
successfully used to image ID correlation functions ob- 
tained at y/s 
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200 GeV [1J|. Briefly, the technique 
numerically inverts the ID Koonin-Pratt equation, 

Ciq) - 1 = Riq) ^AnJ drr^Ko{q, r)S{r) (1) 

which relates the two-particle angle-averaged ID correla- 
tion function, C(g), to the ID source function or image, 
S(r). The latter gives the probability of emitting a pair of 
particles with a separation distance r in the PCMS. The 
ID kernel KQ{q,r) incorporates the effects of Coulomb 
interaction and of Bose-Einstein symmetrization. 

Contamination by uncorrelated pairs (weak decay 
products accepted by the track selection cuts, misiden- 
tified particles, etc.) dilute the correlation and reduce 
R{q). It has been confirmed by simulation that the con- 
tamination is approximately constant in q, so that the 
reduction factor can be assumed to be g-independent. 
The source function S{r) then gets reduced by the same 
r-independent factor due to the linearity of Eq. 

Figure [IJa) shows data points for the ID correlation 
function in relation to the imaged source function in 
Figs. [Ijb) and (c), for mid-rapidity, low pr pion pairs. 
The source function indicates a tail for r > 15 fm which 

is qualitatively similar to that reported for RHIC data 
in Ref.[l3l. a check, the extracted source function is 
used as input to Eq. ([1]) to obtain a restored correlation 
function also shown in Fig. [T](a) ; excellent consistency is 
observed. 

In parallel to the imaging procedure, two different 
functional forms were used to fit the measured correla- 
tion function directly, as discussed below. The conclu- 
sion from the fits (see Fig. [T]) is that a triaxial Gaussian, 
frequently termed ellipsoid, as used in traditional HBT 
methodology, poorly describes the correlation function at 
low (? > 13 MeV/c (Fig. [UJa)), and this leads to a devia- 
tion from the tail of the imaged source function for large 
r > 15 fm (Fig. [Hb)). Fig. [^c) highlights the fact that 

the tail for r > 15 fm contains a very significant fraction 

of the source. On the other hand, the Hump function 
(cf. Eq. ([7]) and discussion below) gives a good fit over a 
more extensive range. 

For systematic access to the 3D source function S{r), 
the 3D correlation function C(q) and source function 
S{r) were both expanded in a series with correlation 
moments R''a^...ai['i) and source moments 5'^^ ..^^(r) in 
a Cartesian harmonic basis representation: 

C(q) - 1 = i?(q) = 5] RU..aMX...a^^^)' 
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FIG. 1: (color on line) Angle-averaged correlation function 
(top panel), source function (middle) and radial probability 
density (bottom) for mid-rapidity pion pairs. Filled circles 
show correlation from direct averaging of the data. Error 
bars indicate statistical errors only; systematic uncertainties 
are smaller than statistical ones. Open circles represent cor- 
relation from fitting the data using angular decomposition. 
Squares show the imaged source and correlation correspond- 
ing to the imaged source. The dotted and solid lines rep- 
resent, respectively, the fitted Gaussian and Hump Eq. (O 
sources and their corresponding correlation functions. 



where / = 0,1,2,..., at = x,y or z, ^qi...^, (f^q) are 
Cartesian harmonic basis elements (fiq is the solid an- 
gle in q space) and Ra^,,,ai{<l) are Cartesian correlation 
moments given by 



(21 + 1) 



-/^<....,(^^.)i?(q) 



(4) 



I Ql-.-Qi 



Here, the coordinate axes are oriented so that z (long) is 
parallel to the beam direction, x (out) points in the di- 
rection of the total momentum of the pair in the LCMS 
frame and y (side) is chosen to form a right-handed co- 
ordinate system with x and z. 

The correlation moments, for each order can be cal- 
culated from the measured 3D correlation function using 
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is similar to the ID Koonin-Pratt equation: 
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FIG. 2: (color online) Correlation moments for multipolarity 
I = 2 (left panels), and I — 4 (right panels) for mid-rapidity 
■k'^-k^ and 7r~7r~ pairs. Error bars indicate statistical er- 
rors only; systematic uncertainties are smaller than statistical 
ones. 



K,...c.M) = 47ry drr^ifz(<z,r)5i^...„,(r), (6) 

but now pertains to moments describing dilFerent ranks 
of angular anisotropy /. Since the mathematical structure 
of Eq. ^ is the same as that of Eq. ([T|), the same imaging 
technique can be used to invert the kernel Ki of the re- 
lation to extract the source moment -S*^^ ^^(r) from the 
corresponding correlation moment B}ai...ai{l)- Finally, 
the total 3D source function is calculated by combining 
the source moments for each I as in Eq. 

Alternatively, the source function can be extracted by 
directly fitting the 3D correlation function with an as- 
sumed 3D shape for the source function. Since the 3D 
correlation function can be represented by the Cartesian 
moments in the harmonic decomposition, the 3D fit cor- 
responds to fitting the six independent non-trivial mo- 
ments simultaneously with a trial source function. 

Figures [l][2] show the result of direct fits to the inde- 
pendent correlation moments with two 3D functions: (a) 
a single triaxial Gaussian, or ellipsoid, (dotted curve) and 
(b) a Hump shape (solid curve). As mentioned, the el- 
lipsoidal fit, with four free parameters, fails to capture 
the low q behavior in C(g) and the large r behavior in 
S(r). On the other hand, the Hump function, with six 
free parameters, gives a good fit. The form of the Hump 
function is 



Eq. (HI). Alternatively, Eq. ^ can be truncated so as 
to include all non-vanishing moments and expressed in 
terms of independent moments only. As expected from 
symmetry considerations, moments odd in any coordi- 
nate were found to be consistent with zero within statis- 
tical uncertainty. Up to order I = 4, there are 6 inde- 
pendent moments: i?^2; R'y2^ ^xA: ^y4 ^Ild Rt2y2^ 

where i?^2 is shorthand for R^^ etc. The independent 
moments can then be extracted as a function of q by fit- 
ting the truncated series to the experimental 3D correla- 
tion function with the moments as the parameters of the 
fit. The present analysis emphasizes the second method, 
with the moments computed up to order I — 4 (higher 
order moments are found to be negligible). The moments 
are shown in Fig. [2l for the multipolarity orders of I =2 
and 4, and in Fig. [TJi for ^ = (1 + R'^iq) = C°{q)) 
The magnitude of the moment R^2y2 comparable to 
that of i?^4 and i?^4 . 

Substitution of C(q) and 5'(r), from Eq. ^ and 
Eq. ([31), into the 3D form of the Koonin-Pratt equation 



C(q)-1 = J drK{q,r)S{r) 



(5) 



results [13] in a relationship between corresponding corre- 
lation i?^^ a; (?) ^i^'i source moments 5*^^ which 



S{x, y, z) — Aexp 



I 



(7) 

where — + y'^ + and the coefficients /s and 
// of the short and long-ranged components are given 
by !s = + (r/ro)^)] and fi = I - f, respectively. 
Here, the argument of the exponential shifts the behav- 
ior from that of a simple spherically symmetric Gaussian 
for r ^ To to that of a triaxial Gaussian for r ^ rg . The 
parameter A regulates the fraction of pion pairs of which 
correlations are described in terms of the Hump function 
(for fit parameter values see Ref. [2^). 

Source imaging involves no assumptions on the ana- 
lytical shape of the 3D source function. On the other 
hand, the moment fitting explicitly invokes a particular 
form for the 3D source function. The ellipsoid fit pro- 
duces a /ndf value of 6.8 while the Hump produces 1.2, 
which indicates a better fit to the observed correlation 
moments, as is visually evident in Figs. [T{a) and[2| Close 
agreement between the experimental data, the Hump fit 
and the restored correlation moments from imaging (see 
Figs.[T][2]) strongly suggests that this assumed functional 
form properly represents the emission source. However, 
the uniqueness of the source function is, for example, not 
guaranteed beyond the region to which data are sensitive 
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such as r > 40 fm or where the source function is very 
small. 

Figures [3t^d)-(f) show comparisons between two-pion 
source functions obtained via the fitting (lines) and the 
imaging (squares) techniques. The ellipsoid fit function 
(dotted line) underestimates the source image (squares) 
and Hump fit function (solid line) for r > 15 fm in the x 
and z directions while the Hump fit function is in good 
agreement with the source image in the cc, y and z direc- 
tions. This consistency check emphasizes the high degree 
of integrity with which the 3D source function is being 
extracted. The source function in the z direction is char- 
acterized by a long tail which extends beyond 30 fm. 
The source function in x also has a non-Gaussian tail, 
which, for this low px cut is less prominent than that in 
z. These aspects are decidedly different from those of a 
RHIC study [H. 

The difference between the source functions from the 
ellipsoid fit and imaging procedures is also evident from 
a comparison of the corresponding correlation functions 
in the x, y and z directions as shown in Fig. [3]Ja)-(c) re- 
spectively. Again there is consistency between the data. 
Hump fit and restored correlation functions in all three 
directions while the differences between the ellipsoid and 
Hump fit sources for r > 15 fm arc manifest by dif- 
ferences between the respective correlation functions for 
g > 15 MeV/c. 

The event simulation code THERMINATOR allows for 
tests of the emission dynamics and of the breakup time of 
the reaction systems jll, [l^ El] . The code simulates 
thermal emissions from a cylinder with input transverse 
radius Pmax- Bjorken longitudinal boost invariance is as- 
sumed, and an expansion with transverse radial velocity 
Vr{p) = (p/Pmax)/(/3/Pinax + «*), whcrc = 1.41, in the 
Blast- Wave mode of the code. A fluid element ring, de- 
fined by p and z, breaks up at proper time r and lab 



The freezeout hypersur- 



frame time t where = r' 
face is specified hy t = tq + ap where a, the space-time 
correlation parameter, is set to -0.5 as was found in Ref. 
[20| . The negative value of a implies "outside-in" burning 
of the source i.e outer particles are emitted earlier than 
inner ones, while a positive value of a would imply the 
reverse i.e source emission from inside out. An emission 
duration parameter At is also needed to achieve a good 
fit. All known hadronic resonance decays are included. 

THERMINATOR parameters Vt T, hb, Us, Mi and 
a are taken from Ref. jH, [l^ 0, US] as obtained from 
spectra and particle yields. Values of Pmax, To and At 
were obtained by matching THERMINATOR's gener- 
ated source function to data shown in Figs.[3fd - f). The 
value of the transverse radius Pmax is chosen so as to 
reproduce the source function profile in the y direction; 
S(ry) is insensitive to tq and At. The proper lifetime tq 
is determined by the short-range behavior of the source 
function profiles in the x and z directions. The proper 
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FIG. 3; (color online) Correlation C{qi) (left panels) and 
source S{ri) (right panels) function profiles for 7r^7r+ and 
7r~7r~ pairs in the outward x (top panels), sideward y (middle) 
and longitudinal z (bottom) directions. The use of symbols is 
analogous to that in Fig. [1] Error bars indicate statistical er- 
rors only; systematic uncertainties are smaller than statistical 
ones. Here, I — 4 moments make negligible contributions. 



emission duration is then determined by the tails of the 
source profiles in the x and z directions. 

The calculation gives a good match to the experimen- 
tal source function in the x, y and z directions with a 
transverse dimension pmax = 7.5±0.1 fm, proper lifetime 
To (t = To at p = 0) of 7.3 ±0.1 fm/c, a proper emis- 
sion duration At ~ 3.7 ± 0.1 fm/c and a — —0.5 (solid 
circles) [28|. The errors quoted are from the matching 
procedure alone. With these values of Pmax, tq and At 
we have reexamined the role of a = —0.5 i.e. outside-in 
burning. Figure 2] shows a comparison of the THERMI- 
NATOR source function, calculated using various values 
of a and other previously tuned parameters [27| , with the 
extracted source function. The open symbols show that 
the calculations with a > overstate the extracted source 
function profile in the z direction. Attempts to com- 
pensate for this overshoot via different combinations of 
Pmax, Tq and At wcrc unsuccessful. Therefore, this fail- 
ure suggests that a negative value for a, hence "outside- 
in" particle emission, is required to reproduce the ex- 
tracted source function. The success of the THERMI- 
NATOR model simulation in precisely reproducing the 
experimental source function indicates consistency with 
approximate boost invariance at mid-rapidity, blast-wave 
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FIG. 4: (color online) Source function profiles, S(n), compari- 
son in the (a) x, (b) y and (c) z directions between the imaged 
data (squares) and THERMINATOR Blast- Wave model with 
various values of a (circles and triangles). Error bars indicate 
statistical errors only; systematic uncertainties are smaller 
than statistical ones. 

dynamics for transverse flow, and outside-in burning in 
the evolution of the expanding system. 

Results from this study and those from Ref.Q de- 
pend on the different analysis techniques and models em- 
ployed. The deduced time scales are similar but the ge- 
ometric transverse radius is quite different. This differ- 
ence results from the inclusion of resonances in THER- 
MfNATOR, as well as different parametrizations of T 
and Vr{p). Conclusions from these THERMfNATOR pa- 
rameters are, of course, model dependent and therefore 
not necessarily unique. Different model assumptions may 
possibly lead to different pictures of the reaction dynam- 
ics m. 

In summary, we have presented a three-dimensional 
femtoscopic study of the two-pion source function in 
Pb-I-Pb collisions at Vs^^y = 17.3 GeV. A model- 
independent imaging/fitting technique reveals prominent 
non-Gaussian tails in the outward and longitudinal di- 
rections of the extracted source function. THERMI- 
NATOR Blast- Wave model calculations, incorporating 
Bjorken longitudinal flow, give a near-exponential tail 
in the longitudinal direction consistent with observa- 
tion. The space-time correlation parametrization sug- 



gests outside-in burning and provides values of the proper 
time for breakup and the emission duration. 
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